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ABSTRACT 

The predicted emission spectrum of N ii is compared with observations of permitted 
Hnes in the Orion nebula. Conventional nebular models show that the intensities of 
the more intense lines can be explained by fluorescence of starlight absorption with 
a N abundance that is consistent with forbidden lines. Lines excited mostly by re- 
combination on the other hand predict high N abundances. The effects of stellar and 
nebular parameters and of the atomic data on the predicted intensities are examined. 
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1 INTRODUCTION 

■ Recent deep spectroscopic surveys of the Orion nebula (Es- 
\ teban et al. 1998; Baldwin et al. 2000, hereafter BVV; Es- 

teban et al. 2004, hereafter EPG) show that N"'"^-|-e recom- 
\ bination cannot explain the intensities of the N ll permitted 
lines with a nitrogen abunda nce that is con sistent with the 
N II forbidden line intensities. ISeatonI |^^^ first suggested 
, the possibility that permitted lines of C and O ions in neb- 

■ ulae may be excite d by continuum fluorescence of starlight, 
' and lOrandil (|l97(]t ) noted that this also must be an impor- 
, tant mechanism to excite the Nil lines in Orion. Grandi 
' suggested that additional absorption of photons of the He I 
; Is^ ^So-ls8p ^P? line at A 508.643 A from the diffuse field 
, of the nebulae by the N ll 2p^ ^Po-2p4s ^P? transition at 

' A 508.668 A followed by decay to 3p terms would enhance 
the observed intensity of some lines. There i s few d irec t evi- 
dence of the plausibility of this mechanism. IBVVI and lEPGl 
observed lines at A 3829.92, A 3838.47 and A 3856.27 that 
could be produced by the 3p ^Pi-4s ^Pj, 3p ^P2-4s ^Pj 
and 3p '^P2-4s transitions respectively. Unfortunately 
those lines are probably blended with lines from other ele- 
ments, and their ide ntification and intensity are uncertain. 
ISharpee et al.l i200 have detected those and other lines of 
the 3p ^P-4s ^P° muhiplet in a planetary nebula (IC 418), 
where the N ll spectrum is probably excited by fiuorescence. 
The 4s ^P§ level requires pumping of the the Nil 2p P2- 
2p4s ^PjA 508.697 A transition, which lies 32kms"^ from 
the He I line. Other lines from the 4s ^P° term like those 
of the 3p ^D-4s ^P° muhiplet in the AA 3311.42-3331.31 
interval, which should be as intense as the 3p '^P-4s ^P° 
multiplet, were not detected by EPG. The efficiency of this 
Bowen-type line fluorescence depends heavily on uncertain 
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nebular parameters that are needed in the theory of line ra- 
diative transfer, and modeling can become quite arbitrary. 
Therefore we will not consider it in this wor k (for further 
discussion see lEscalantdliooi iLiu et aljEool . 

The critical parameters that determine the intensities 
of the lines and the relative importance of the fluorescence 
mechanism over the recombination process are the and 
N"*"^ column densities in the gas and the stellar UV radia- 
tion field. Absorption of a UV photon by a resonant tran- 
sition between a ground configuration state and an excited 
state has a higher probability of subsequent reemission in 
the same transition. Decay to an intermediate state will be 
favoured when the optical depth of the resonant transition 
is large, and the resonant photon is scattered a few times 
until it is converted into a lower energy photon producing 
a subordinate line. However the optical depth of the reso- 
nant transition must not become too large in order to allow 
enough resonant photons to penetrate into the N+ zone. 
The efficiency of the continuum fiuorescence excitation in 
N""" depends on the transfer of resonant photons between 
the lowest resonant transition that can produce a subordi- 
nate line, 2p^ ^P-3d ^D°AA 533.51-533.88 A, and the ioniza- 
tion limit at 419 A. We have also included transitions from 
the 2p'^ ^^02 and ^So metastable states-populated mostly by 
coUisions-to other singlets in order to consider the observa- 
tion of singlet lines in Orion. In the singlet system the lowest 
transition that produces a subordinate line is 2s^2p^ ^Sq- 
2s2p^ ^P?A 745.84 A. Fig. shows the observed transitions 
in Orion of the singlet and triplet N 11 systems. Possible ob- 
servations of quintet lines and higher excitation states are 
discussed in section ITU 

The observed intensities of permitted lines in planetary 
nebulae are often higher than their intensities predicted by 
recombination r ates with CNO abundances measured from 
forbidden lines llLiu et al.lll995l. 1200 ll : iLuo et aljl200ll and 
references therein). Some of the Nil lines observed in PNe 
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Figure 1. The singlet and triplet Nil systems. Terms 3d ^P°, 
3d 3d° and 4s ^P° are connected by resonant transitions to the 
2p2 3p ground term at eV. Wavelengths are term— averaged 
values in A. 



have 4f upper levels, which are excited mainly by recom- 
bination. In Orion there exists a similar situation and we 
will show that fluorescence cannot account for the excess 
intensity of lines from 4f levels. The accuracy of the recom- 
bination theory can be more easily tested by comparing line 
ratios from 4f decays because all the transitions involved are 
optically thin and the absolute line intensities have a similar 
dependence on temperature and density. 

The fluorescence theory is more difficult to test because 
it depends, often non-linear ly, on several model parameters. 
We will use the ratio of predicted over observed intensi- 
ties, averaged over all N ll lines with reasonably accurate 
identifications and measurements, hereafter referred to as 
R = {Icai/Ioba), to test the model parameters. The scatter 
around 7? will be used as a test of the atomic data and the 
details of the transferred stellar continuum. 

Realistic nebular models and hot star model atmo- 
spheres along with available atomic data bases can ex- 
plain successfully the intensities of a majority of the for- 
bidden lines in Orion as well as general observed trends 
of ion abundances in Galactic HII regions (see for example 
iBaldwin et alll99ltlstasiiiska fc Schaereill997ft . This paper 
demonstrates that the intensity of most of the N ll permitted 
lines in the Orion nebula can be predicted by fluorescence 
of the starlight continuum and some contribution of recom- 
bination by these models with currently accepted physical 
conditions and abundances of the nebula. 



2 CALCULATION OF POPULATION 
DENSITIES 

2.1 Atomic processes 

The important processes that can produce excited states 
with low principal quantum numbers in N"*" at nebular 
temperatures are absorption of UV photons by transitions 
from ground and metastable states and recombinations of 
N"*"^. Subsequent decays of these states produce the optical 
Nil spectrum. UV radiation is provided mostly by the star 
C Ori. The contribution of the diffuse continuum to the 
absorption is negligible. 

The number density of a excited state j at a certain 
point in the nebula in steady state, rij, is given by 



3 



kj 



(1) 



k>j 



where Aj = Aji is the spontaneous decay rate, aj is 
the recombination (radiative plus dielectronic) coefficient to 
state j, Ug is the population density of one of the five N"'" 
ground and metastable states, 2p^ ^Po,i,2,^ D2,^ So, and Pgj 
is the radiation absorption rate. If is the local starlight 
mean intensity averaged over the absorption profile. 



f3g 



hv 



l9i 
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tion and fgj is the f-value of transition g ^ j- Jv is atten- 
uated by continuum opacity and geometry. The "pumping 
probability" to account fo r the att enuation due to the tran- 
sition 3 — > j as defined bv lFerlandt (.1992.') is 



7gj 



-To<t>(x) 



(j)(x) dx 



(3) 



where 4>{^) is the normalized Voigt profile for a displacement 
X = {v — i/q)/ l^v, and tq = OgjNg/ l^v is the mean optical 
depth for a column density Ng — J Ug dr integrated along a 
ray from the star. 

An actual steUar spectrum shows absorption lines su- 
perposed on the continuum and the P-Cygni profiles of the 
wind. If the structure of the stellar absorption spectrum 
varies in scales comparable to the width of AA ^ 0.01 A 
of the UV resonant lines in the nebular gas, equation ((31 
should be changed to 



lai ~ I '<P{x)e ^°'^'^'(j!)(a;) da; 



(4) 



where tp{x) = Ji,/ Jv is the continuum stellar profile around 
the resonant line. Grids with that resolutio n are avail- 
able in the optica l jMurphv fc MeiksinI l2004ll In the UV 
ISmith et all (|22q3) and ISternberg et alJ (|20o3^ have pub- 
lished models with a lower resolution. 

The populations of the metastable states at nebular 
temperatures are controlled by electron collisions, and are 
nearly independent of other excited states. The system of 
equations Q is triangul ar, and can be solved in terms of 
the cascade matrix, Ckj llSeatorJll959t) . Ckj is the probabil- 
ity that a state j is produced by the excitation of state k 
followed by radiative decays by all possible routes ending in 
j. Equation Q thus becomes 
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where 

a'f — atCkj 



k>j 



(5) 

(6) 
(7) 



Table 1. / absorption values for some Nil resonant multiplets. 



are the effective recombination coefficient and the effective 
fluorescence rate respectively. In practice the sums in equa- 
tions Q and Q must be truncated and a correction must 
be added to equation ^ by extrapolation or by using hy- 
drogenic populations (Escalante & Victor 1990, hereafter 
EV) because recombination coefficients of individual levels 
decrease slowly with n and the contribution of states with 
high angular momentum and high n must be included. Care 
must also be taken to avoid roundoff errors in the summa- 
tions. The effective fluorescence rate in equation Q is less 
sensitive to that type of correction because the absorption 
rate Pgk is non-zero only for states connected to the ground 
and metastable states by dipole-allowed transitions, and 
decreases more rapidly with the principal quantum num- 
ber n of the upper state. The contribution of levels with 
4 < n < 12 and orbital angular momentum < i < 2 in 
equation is less than 5% for the N ll lines observed in 
Orion. The contribution of f states to the fluorescence rate 
is even less important because they are not connected by 
resonant transitions to the ground term, and the transitions 
nd mf have a low relative probability. By eliminating 
states with n > 9 and i > 2 in equation Q, the CPU time 
decreases by a factor of 10. However the f states (and higher 
I states) must be included in the cascade due to recombina- 
tions, which tend to favor high-angular momentum states. 

Most of the observed N ll transitions in Orion come from 
decays of the 3p and 3d triplet terms. Practically all the ex- 
citations of the 3d terms are produced by direct absorptions 



^P-2p3d ^P°AA 529.36-529.87 



in the multiplets 2s^2p^ 
and 2s^2p2 3p-2p3d 3D°AA 533.51-533.88 A. Cascades con- 
tribute negligibly to the populations of the 3d terms. The 
3p terms, not being connected by direct transitions to the 
ground term, have more varied excitation channels. Between 
half and 80% of excitations of 3p terms come from decays 
of 3d terms. The rest comes mostly from absorptions in 
the multiplet 2s22p2 3p_2p4s ^P" A A 508.48-509.01 A and 
higher s and d states. Table shows a comparison of /- 
values for these multiplets, which are critical in our calcula- 
tions. 

The advantage of using the cascade matrix is that it 
needs to be computed once in the model in either the op- 
tically thin case (case A) or the extreme thick case (case 
B) because it depends solely on the Einstein coefficients, or 
more precisely, on the branching ratios Pkj — Akj / ^ . Aki- 
In the escape probability formalism the equations must be 
modified by substituting the Einstein coefficients of tran- 
sitions connected to the states of the ground configuration 
2p^ ^P, ^D, by P^Akg, where is the escape probability 
to be discussed below. Pe is a local quantity that depends 
on the optical depth of the line as a function of position and 
thus the cascade matrix must be recomputed at each point 
in the nebula. However only a small fraction of the matrix 
elements depend on Pe. We found that recomputing only 



Multiplet 


A (A) 
(1) 


(2) 


/ 

(3) 


(4) 


(5) 


2p2 3p-4s 3p° 
2p2 3p-3d 3p° 
2p2 »P-3d 


508.74 
529.68 
533.67 


0.00987 

0.102 

0.294 


0.0175 

0.175 

0.575 


0.0105 

0.103 

0.300 


0.0124 

0.138 

0.378 



(1) 

(2) 
(3) 
(4) 
(5) 



Multiplet-averaged wavelengtl 



Wiese et all il99i 



Victor fc Escalantd J198I 



Seaton_et_alJ (199^ (O pacity Project). 



NussbaumCTfeStoetl 11984) . 



this fraction of the cascade matrix reduces the CPU time by 
a factor of 30. 



2.2 Escape probabilities 

We have used the escape probability theory to handle nearly 
200 resonant transitions that produce 2630 transitions by 
fluorescence. A review of the limitations of this theory has 
been given by iDumont et al.l i2003t) , but a full line trans- 
fer calculation is beyond present computation capabilities. 
There are many approximations for the escape probabil- 
ity function, which differ by several factors at high optical 
depths. In an ionized nebula absorption of reso nant pho- 
tons by an overlapping continuum is important iHummej 
Il968h . In a uniform slab of total optical thickness T the es- 
cape probability at depth r is Pe = bF{b) + {1/2)[K2{t, b) + 
K2(T — r, b)], where the functions F and K2 are defined in 
iHummer fc Storevl il992h . and b = kc/ki is the continuum- 
to-mean line opacity ratio. The term 



bF{b) = 



ki + k, 



■4>{x) dx 



(8) 



is the probability that the resonant photon will be lost by 
continuum absorption. The effect of Pe is to increase con- 
siderably the probability that a resonant absorption decays 
into the subordinate line rather than reemitting the reso- 
nant photon. The average ratio, R — (/cai//obs), increases 
by a factor of 20 between the two limits Pe = 1 (case A) 
and Pe = (case B) respectively. A spectrum dominated by 
fiuorescence, however, must be in an intermediate regime in 
order to allow the penetration of UV stellar photons through 
a large column density of absorbers at the same time that 
the resonant photons are scattered repeatedly. 

All the resonant photons that produce the Nil optical 
spectrum have a high probability of being lost by conversion 
into other lines besides being absorbed by the continuum, 
and the probability of a large number of scatterings is small. 
Consequently we used a Doppler core in the calculation of 
K2{r,b) and F{b). 



2.3 Atomic Data 

The main sou rce of A- an d /-va lues for this work is the 
compilation of I Wiese et al.l lll996l) , which is based pr imar- 
ily on the Opacity Project data dSeaton et aljll994f) and 
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Table 2. Effective fluorescence rates for some Nil lines. 



Line 



MA) P..E,/3^« 



(10- 





i 


- j 




(1) 


(2) 


(3) 


(4) 


3s 


3po 


-3p 


4630.5 


10.8 


25.7 


16.1 


13.7 


3p 


^D3 


-3d 


4803.3 


4.6 


4.5 


4.5 


4.6 


3s 


3po 
^2 


-3p ^Da 


5679.6 


7.1 


16.3 


10.3 


8.9 


3p 


3P2- 


-3d 3d° 


5941.7 


8.0 


8.6 


8.6 


8.0 



Wiese et all 



Victor fc Escalant j <198 



(1) Only 

(2) Only 

(3) Same as (2) plus doubly-excited configurations 
and AS 7^ transitions. 

(4) All transitions included. 



tween the three data sets is in the 4f terms, where most 
of the contribution to the recombination comes from levels 
with small non-hydrogenic effects. The main uncertainty is 
in the line fractions involving 4 f ter ms, where LS-coupling 
is not a good approximation. IEVI used LK coupling for 
these terms, and line fractions fo r other couplings are avail- 
able llEscalante fc Gongorall99oll , but general intermediate- 
coupling calculations for those states are clearly needed. 

This work uses the recombination coefficients of 
I EVI with branching rat ios given by the A- values of 
IVictor fc Escalantd Il988l) and lwiese et alj Jl99d) . Effective 
recombination coefficients for the levels of a term were ob- 
tained by assuming that the coefficients are proportional to 
the statistical weights of the levels. 



configuration-interaction calculations, but also includes in- 
termediate coupling calculations and experimental measure- 
ments. That compilation has data for most lines of the se- 
ries 2s^2pCP)nl with n < 5 and Z < 2, and 2s2p^(*P)nZ 
with n < 2 and / < 1, including spin-forbidden transitions 
between quint et and triplet terms. The m odel potential cal- 
culations from lVictor fc Escalantel (ll98Sh were used for the 
rest of the transitions needed in the cascade matrix of the 
effective fluorescence rate in equation Transition prob- 
abilities between fine structure levels wer e obtained by ap- 
plying LS fractions to the multiplet data llAllerJll973ll . 

Table|5|shows the effect of different atomic databases on 
the effective fluorescence rate Q summed over the ground 
and metastable states times the branching ratio of the sub- 
ordinate line. The continuum is a blackbody spectrum at 
T = 37kK with a photon flux of 9.24 x 10"" cm"^ s"^ Hz"^ 
at A 533.7 A and ro = for all lines, i.e., 7o, = 1, 
Pe = 1. The compilation of IWiese et alJ lll996ll does not 
include transitions to levels with principal quantum num- 
ber greater than 5 and produce lower fluorescence rates of 
fines with upper p levels because those levels have cascade 
co ntributions from higher lev els. The model potential data 
of I Victor fc Escalantd (|^8^ does not have transitions to 
doubly-excited conflgurations and consequently give higher 
rates as shown in the fourth column of table |21 Transi- 
tions between states with the 2s2p^ configuration and the 
singly excited configurations 2s^2p3p can change the branch- 
ing ratios significantly. The entries in the fifth column have 
been complemented with transitions to states with the 2s2p^ 
core configuration a nd spin-forbidden transitions taken from 
I Wiese et al.l l)l99fi|) . The last column combines both data 
sets. The cascade matrix elements tend to be similar for dif- 
ferent data sets because systematic differences in the atomic 
parameters between different data sets cancel out in the 
branching ratios Pji = Aji/Aj. 

A recent calculation of effe c tive r ecombination co- 
efficients by iKisielius fc StorevI (|2Qq3) shows a gen- 
eral agreement with the model potential calculations of 
lEscalant e fc Victor] l|l99 (t) (he reafter EV) and the calcula- 
tions of ^iauimo^^Lri|l99ir i for the 2s22p(2p°)3d and 4f 
configurations. The most important differences between the 
three calculations are in the branching ratios of multiplets 
from the 2s^2p('^P°)3p '^D term. The accuracy of the model 
potential in this case was limited by the lack of observed 
energies in the 2s^2p(^P°)np series. The best agreement be- 



3 MODEL CALCULATIONS 

3.1 Nebular models 

In order to determine the electron, N"*"^, and densities 
in equation 0, as well as the temperature and opacity 
at each point in the nebula, we used th e codes cloudy 
jFerlandlll99d. version 9 0.05) and nebu jPeauignot et all 
l200ll:lMorisset et al.l2002ll . Models of the Orion nebula sup- 
port the existence of a main emitting layer at the back of 
a cavity in the OMC-1 molecular cloud. T he thickness of 
the la yer is highly variable across the nebula ijWen fc O'DelJ 
Il995h . We approximated the layer by a plane-parallel model 
at constant pressure with CLOUDY and constant density with 
NE BU. CLOUDY allowed US to use the grain composition used 
bv lBaldwin et all il99ll) in their Orion model. Predicted line 
intensities by CLOUDY show some sensitivity when the radi- 
ation is included in the pressure law (Baldwin ct al. 199(|). 
Models with a constant gas pressure produce larger col- 
umn densities than models with a constant gas plus radia- 
tion pressure. We have not tried to find a single model fit to 
the observed forbidden line intensities in Orion. Instead we 
have run a series of models to find the most important de- 
pendencies of the N II lines excited by fiuorescence on model 
parameters. 

3.2 Model atmospheres 

The Orion nebula is mostly excited by C Ori. The other 
Trapezium stars increase the fiuorescence of the N 11 lines by 
less than 2% and will not be considered in this calculation. 
6^ C Ori is a class V star with variable wind features that 
produce uncertainties in the determination of its spectral 
classification and effective temperature. Different authors 
give values close to Teg = 39 kK and log(o) = 4 for this 
star jHowarth fc Prinialll989l: iHiflenbrandl 119971 and ref- 
erences therein). However comparisons of the intensities of 
forbidden lines with nebular model predictions suggest tem- 
pe ratures as low as 36 kK (BVV). Metalicity measurements 
bv lCunha fc Lambert! (Il994l) show that the Trapezium stars 
are slightly underabundant with respect to the Sun. 

The calculation of the fiuorescence of the gas needs a 
high resolution stellar spectrum. This is important with hot 
massive stars, which have expanding atmospheres with a 
dense forest of absorption lines and broad overlapping P 
Cygni profiles. The emission peaks and absorption troughs of 
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Table 3. Model parameters for the basic model (BM) 





37 kK 


n 


10" cm- 


~3 


log 3 


4.0 


He 


0.095 




Z 


IZq 


C 


3 X 10- 


-4 




1012.9 s- 


In* 


6 X 10" 


-5 


e 


0° 


o 


4 X 10" 


-4 


turbulent 


8kms-i 


Ne 


6 X 10- 


-5 



* Variable parameters in this work. 



the overlapping P Cygni profiles can increase or decrease tlie 
absorption rate in equation ^ by large factors and ciiange 
tiie fluorescence excitation rate significantly. On the other 
hand nebular models usually smooth the spectrum of the 
exciting star to calculate the ionization structure. In order to 
take into account the detailed spectral structure of the model 
atmosphere, we input the unattenuated stellar spectrum into 
the nebular model and extracted the predicted attenuated 
local continuum at each point in the nebula. To calculate the 
pumping rate in equation Q, the full resolution spectrum 
was read and interpolated at the absorption frequency and 
was scaled by the ratio of the attenuated to the unattenuated 
stellar continuum predicted by the nebular model. 

The motion of the star with respect to the gas intro- 
duces a Doppler shift that can change the intensity of the 
continuum at the absorbing frequen cies. The proper mo - 
tion velocity of 6^ C Ori is uncertain llWen fc 0'DelJll995l) . 
Doppler shifts of up to ±30kms~^ did not change predicted 
line intensities by more than a few percent with the reso- 
lution of about AA ~ 1 A in the far UV of the model at- 
mospheres that we used. Therefore we assumed a static star 
with respect to the gas. 

Recent model atmospheres of O stars include the effects 
of line blanketing and line blocking of the stellar wind. We 
used model atmospheres calculated with the wmbasic code 
jPauldrach et alJl200 it [Sternberg et"alll2003ri to account for 
these effe cts. We also us ed the LTE, line-blanketed atmo- 
spheres of lKuruc3 il99ll) for comparison purposes. 



4 BASIC MODEL 

All transitions of the observed permitted lines in Orion end 
in excited states and are optically thin. Their intensities can 
thus be obtained by integration of the emissivity along the 
line of sight: 

J ^J""^ '^'^ (9) 

We now examine the dependence of the fiuorescence excita- 
tion on the stellar spectrum and th e de nsity and compare 
them with the observations of iBVVl and lEPGi 

We have adopted a set of central values for the param- 
eters of the models that are given in table |31 Those values 
are close to the ones recommended by Baldwin ct al. ( 1991,), 
and we will refer to them as the basic model (BM). We dis- 
cuss below a few features of this model. 



4.1 Gas abundances and density. 

The abundance in the nebula can be estimated from 
the measurements of lines with upper 4f levels, which are 
populated mostly by recombination. Table |1| lists the ob- 
served intensities of these lines in Orion and their pre- 
dicted recombination emission rates normalized to the ob- 
served intensities and recom bin ation rate of the 3s ''P2- 
3p ^DaA 5679.56 A line. lBVVl and lEPGl measured the lines at 
AA 4236.91, 4237.05 and 4241.784, which account for 92% of 
the total strength of the 3d ^D°-4f F muhiplet if the 4f lev- 
els are described by an LK coupling scheme as shown in the 
table. We do not include in the abundance estimation the 
A 4242.49 line because its measured intensity is much higher 
than the one predicted from the recombination theory. At 
Te = 10* K the effective recombination coeff icient of multi- 
plet 3d ^D°-4f F is 5.6 x 10"^* cm^ s"^ SeW^) . which implies 
N+VH+ = 9.5x10"^ if the lines were produced solely by re- 
combination. This abundance is higher than the one implied 
by the forbidden lines ft; 6 x 10~^ in Orion and many galactic 
Hll regions (e.g- JShaver et al.llT983 : Afllcrbach ct al. 197^. 

We notice that EPGf detected one line from the 3d ^F°- 
4f G multiplet (A 4041.31) although other lines of this mul- 
tiplet theoretically should be more intense than the ob- 
served lines of the 3d ^D°-4f F multiplet. The corr esponding 
N+^/H "*" abundance from this line is 8.5 x 10~^. iLiu et alJ 
J200ltl has reported a abundance of 4.47 x lO"'^ in Orion 
from the AA 4041.31 and 4043.53 lines, which is more con- 
sistent with the abundance from the forbidden lines. 

Several authors jLiu et alJ l2000l . l200lt iTsamis et all 
l2003l:IPeimbert et alj|2004) have detected Nil lines in plan- 
etary nebula with relative rates that are more compatible 
with the recombination theory and will be the subject of 
future work. In a typical H 11 region there is not enough N^^ 
column density to produce a N""" recombination spectrum. 
Fig-Elshows separately the intensities given by equation Q 
(integrated from right to left) due to recombination and fiuo- 
rescence as a function of depth for a line. While the recombi- 
nation intensity grows linearly with distance the fiuorescence 
intensity grows more rapidly. The fluorescence excitation is 
favoured by the more intense starlight continuum and low 
opacity in the near side of the nebula, while the recombi- 
nation emissivity is more uniformly distributed. In the far 
side of the nebula the high optical depth in a resonant line 
scatters more photons and increases the probability of re- 
absorption in the line, thus increasing the pumping of the 
fluorescence emission. The measurements bv iBVVi show that 
most Nil permitted lines are blueshifted by 2 or 3 kms~^ 
with respect to the N 11 forbidden lines, and support the idea 
that they are formed in different layers of the nebula. 

We h ave chosen as a startin g point the set of gas abun- 
dances of iBaldwin et all il996fl . which were derived to flt 
forbidden line intensities with CLOUDY. The only exception 
was nitrogen for which we took a val ue of 6 x 10~^ as a n av- 
erage from the determinations of Ba ldwin et all |^9^ and 
IEPGI Differences in the abundances of the other elements 
between the two sets are not important in N 11 permitted 
line intensities. 

The N^ column density remains remarkably constant 
with varying gas density for a given effective temperature of 
the star, T^s. For the T^g = 37 kK, logg — 4 atmosphere, 
iV(N+) = 2.5 X 10^^ cm"^ for gas densities between 4,000 
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Table 4. Intensities of lines from 4f levels in Orion and PNe [/(N llA 5679.56) = 1]. 



LK Multiplet 


Line 


A (A) 


% 


Theory 


Orion 




PNe 




3d 3L°-4f L'[K]j/ 


J-[K]J' 


J-[K]J' 


(1) 


(2) 


(3) 


C4) 


(5) 


(6) 


(7) 


3d 3Fj-4f G[K]j/ 


3-[9/2]4 


4026.08 


13.9 


0.24 




* 










2-[7/2]3 


4035.08 


23.8 


0.41 






0.24 


0.23 


0.32 




4-[9/2]5 


4041.31 


40.7 


0.70 




0.30 


0.58 


0.53 


0.58 




3-[7/2]4 


4043.53 


17.4 


0.30 






0.35 


0.21 


0.23 


3d 3D5-4f G[K]j, 


3-[7/2]3 


4201.35 








0.14 








3d F[K]j/ 


l-[5/2]2 


4236.93 


20.0 


0.20 \ 


0.23 


0.16 


( o.n 


0.13 


0.16 


2-[7/2]3 


4237.05 


12.7 


0.13 J 


\0.26 


0.20 


0.23 




2-[5/2]2 


4241.24 


3.7 


0.04 






0.04 


0.04 


0.04 




2- [5/2]3 

3- [7/2]4 


4241.76 
4241.79 


16.9 
42.9 


0.17 \ 
0.42 J 


0.20 


0.28 


0.49 


0.54 


0.51 




3-[7/2]3 


4242.49 


1.6 


0.02 




0.28 









* Blended with He I. 
(1) Line fraction in LK coupling. 
(2) Calculated from recombinati on rate s at T = 10* K in case A, 

(3) 

(5) NGC 6153. iLiu et"an l2m ' 



EPG 



(6) M 1-42, 

(7) M 2-36, 



intensities from scanned spectrum. 



Liu et al 



Liu et al 



2001). 
200A. 




r (pc) 

Figure 2. Ion fractions of N+, N+^, fluorescence and recom- 
bination contribution to the intensity of the 5941.65 A line and 
stellar continuum (long dash) at 533.73 A as functions of ge- 
ometrical depth in parsecs. The intensities are normalized to 
the total intensity of the line, 4nl = 5.14 X 10~* erg cm~^ s~^ 
and to the stellar flux on the illuminated side of the nebula, 
AttJi, = 1.402 X 10-3 erg cm-2 s'^ Hz"!. Earth and 9^ C Ori are 
to the left. 



and 25,000 cm ^. Therefore we will adopt a fixed density of 
n = 10" cm-^ 



4.2 Stellar continuum and geometry 

The blister model, in which 9^ C Ori is near the wall of 
the Q MC-1 molecular cloud llZuckermanlll973 : iBalick et alJ 
Il974) . gives the most likely geometry for the nebula. The 
Lyman photon flux of the star at a distance ro from the il- 
luminated face of the nebula is (f)o = Qo/'^htq, and can be 
constrained by the H/3 observed brightness up to a geometri- 
cal factor that depends on the angle of illumination 6 of the 
slab of gas as described by Wen & O'DcU ( 1995). The value 
(j) — 10^^'^ cm"'^ s"'^ and = Q, flxes t he H/3 intensity a round 
AnI = 2.6ergcm"^s"^ observed bv iBaldwin et al.l (ll99J) 
if corrections are considered to account for reflected optical 
and absorbed UV radiations jFerlandl200ll) . The reff=37kK 
is too low for the probable spectral type of 9^ C Ori, but was 
chosen because it reproduces the absolute intensities of the 
fluorescence lines with respect to the H/3 flux and the N 
abundance measured from forbidden lines as discussed in 
section |5. II 

4.3 The recombination spectrum. 

The 4f levels in N"*" are in a n intermediate coupling be- 
tween LK and jK couplings dCowan fc Andrew! Il965() . In 
both couplings the total angular momentum is J = K±l/2. 
The line fractions for either coupling can be obtained from 
lEscalante fc Gongoral lll990l) . We employed an LK classifi- 
cation to use the recombination rates of lEVl 

Table^Jshows that all but two of the observed ratios of 
intensities of lines from 4f levels with respect to lines from 
3d levels tend to be twice as strong and are closer to their 
predicted values by the recombination theory in PNe than 
in Orion. The two exceptions are the A 4242.49 line, which 
seems much stronger than its predicted recombination rate, 
and the line of the 3d ^'D°-Ai G[7/2] multiplet, which vio- 
lates the AL = 0, ±1 selection rule in L K coupling due to 
stron g mixing with the 4f F[7/2] term (|c owan fc Andrew! 
Il965!l . and has no predicted recombination rate by !evI 
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iLiu et all i2000l . I2OOII) observed lines from other 4f terms 
in PNe, which have not been detected in Orion and are not 
listed in table0] The higher intensity of the lines from 4f lev- 
els with respect to the lines from 3d levels in PNe indicates 
that fluorescence is less important relative to recombination 
in these objects because recombination can efficiently popu- 
late levels with high angular momentum while fluorescence 
is limited to resonant levels. One important exception to this 
trend is PN IC 418, which shows a strong enhancement of 
lines f rom 3p and 3d levels with respect to lines from the 4f 
levels iSharpee et aljEoOSl) . This nebula has a lower ioniza- 
tion level than most PNe that favors the excitation of N ll 
lines by fluorescence over recombination as in Orion, and 
will be the subject of a forthcoming paper. 

Fluorescence can only contribute to population of the 
4f levels through transitions from higher d levels, but this 
process is inefficient because the transition probabilities 
are much smaller than transitions to p levels. For exam- 
ple, fluorescence can populate the 4f F levels through ab- 
sorptions in the multiplet 2p^ ^P-5d ^D°A453, but in a 
pure case B the 5d '^D° term decays to the terms 3p ''P, 
3p ^D, 2s2p^(*P)3s ^P, 4p ^P and 4f F among others with 
branching ratios P(5d, 3p ^P) = 0.70, P(5d, 3p ^D) = 0.14, 
P(5d, 3s) = 0.09, and P(5d, 4f) = 0.04. As shown in sec- 
tion l4.1l N"*"^ abundances obtained from recombination rates 
of the 4f lines in Orion are 1.6 times the N abundance de- 
rived from coUisionally excited lines. Suppose that the in- 
tensity of a 3d-4f line were enhanced with respect to If/3 by 
fluorescence populating the 5d term with rate B{5d): 

— P(A) (10) 

where £;(H/3) = / Qg*^n(H+)n, d^, E ^ J alfpn{N+^)ne d£ 
and P(A ) is the branching ratio of the 3d-4f line with 
A given in A. To account for the excess abundance 
of 0.6 obtained from the effective recombination rate 
of multiplet 3d ^D°-4f FAA 4236.91-4247.20, we need 
B(5d) — 0.6i?/P(5d, 4f), which means a rate of pop- 
ulation of the 3p ^P term of B(5d)P(5d, 3p ^P) = 
10. 5i?. The contribution to the intensity of a multiplet 
like 3s ^P-3p ^PAA 4601.48-4643.08 due to this addi- 
tional excitation with N+^/H+ « 6 x 10"* would be 
S(5d)P(5d,3p 3p)P(3p ^P,3s ^P°)/S(H/3) = 7.7x10"'' 
where we took P(3p ^P, 3s ^P°) = 0.36 and an effective 
recombination coefficie nt of the 4f F term a'^p = 1.0 x 
10"" cm^ s"^ at 10* K iIevI) . The strongest line of the mul- 
tiplet, 3s ^P^-3p ^P2 A 4630.54, has an LS line fraction of 
11.25/27 jAllenlll973h and the corresponding increase in in- 
tensity would be at least 0.034 (/(H/3) = 100), which is com- 
parable the observed intensity of 0.048 (EPG) produced by 
more direct cascade routes following absorption of photons 
at AA530 and 534 A by 3d states. Thus the excitation of 4f 
states by fluorescence would produce lines from 3p and 3d 
states with intensities much higher than the observed values 
unless the stellar continuum had an unusual shape that se- 
lectively excited states above the 4f states. Our calculations 
show negligible contribution of fluorescence to the excitation 
of the 4f levels because the absorption rate is much less for 
higher resonant levels than for the 3d leve ls, and point to 
oth er mec hanisms to excite them iTsamis et alji2004ri . 

lEPGl also observed the lines at A 5001.14 and 5001.48 
with upper levels 3d ■'F2^3. The most intense component 



of the multiplet at A 5005.15 is blended with the [O ill] 
line. As with the lines with upper 4f levels, the A 5005.15 
line is mostly excited by recombination because its upper 
level 3d '^F4 can receive only indirect contributions from 
the fluorescence excitation of higher levels. The other lev- 
els, 3d ''F2,3, are connected to the ground sta te through 
weak dipole-allowed transitions jBefl et al Jl 19951) ■ and have 
a substantial fluorescence contribution. Our model calcula- 
tions show that fluorescence contributes less than 5% to the 
intensity of the lines produced by the 3d '^F4 and 4f levels in 
Orion, and it is not sufficient to explain the discrepancy be- 
tween the abundances determined from recombination and 
coUisionally excited lines. 

4.4 The fluorescence spectrum 

Table|K|gives the predict ed i ntensi ties for the observed upper 
terms in Orion bv lBVVl and lEPGl that are excited mostly by 
fluorescence. Although the two data sets are from different 
parts of the nebula, it is important to notice that both sets 
give similar measurements of the N ll permitted line intensi- 
ties with respect to H/3. Uncertain observed intensities due 
to blends, low S/N or dubious identifications are marked 
with "?" in the table as indicated by those authors. We 
also added a question mark to the line at 4987.4 A, which 
is probably blended with the [Fe III] A 4987.20 line, and thus 
has an observed intensity much higher than our prediction. 
The line at 4994.4 A belonging to the same multiplet should 
be theoretically more intense, contrary to the observations. 

The most intense fluorescence lines are triplets con- 
nected by resonant transitions to the ground term, 2p^ '^P. 
At temperatures characteristic of H ll regions, the fine struc- 
ture populations of the ground term are approximately pro- 
portional to their statistical weight, and consequently the 
relative intensities of all the other triplet levels are given by 
the LS line fractions jAllenlll973ll . 

The lines at 5001.15 and 5001.48 A arising from the 
3d ^F° term are blended. We split the total intensity accord- 
ing to the LS line fractions: 21.0:31.1, in order to compare 
them with our predictions in table |5] 

Multiplet 3p ''P-4s ''P° has many strong lines in the 
A 3838.37-3856.06 A interval that were not detected because 
of the lower S/N at the blue end of the spectrum and because 
they are blended with other lines. Therefore they are not 
liste d in t able |S] 

IEPGI marginally detected the 3d ^P2-4p ^Si A 6809.99 A 
line. Lines from the 4p levels must be excited by cascades 
from upper levels in a fluorescence spectrum. We found very 
low intensities for all the lines from 4p levels. The most 
intense line of this type is 2p^ ^P2-4p ^SiA 1060.2 A with 
an intensity of 4 x 10"* with respect to H/3. Our predicted 
intensities of lines from 4p levels in the optical are below the 
instrumental sensitivity. 

In table |S] we also list some singlets because iBVVl 
marginally detected the 2p^ ^D^-3p ^Pi A 4895. 11 A line. 
Singlets can be excited by absorptions from the '^D and '^S 
terms of the ground configuration, spin-forbidden transi- 
tions from the triplets and recombination. The quintet lines 
3p *D?-3d *PiA 4815.62 A, 3s *Pi,3-3p *Di,4A 5535.36 A 
and 3s *P3-3p ^DaA 5551.99 A measured bv .BVV. and .EPCl 
are not listed. These lines have the 2s2p^ ""P excited core 
configuration, and can only be populated with transitions in 
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Table 5. Predicted intensities of the basic model (BM) and ob- 
servations (/(H/3) = 10*) for lines excited by fluorescence. The 
nebular model was calculated with CLOUDY. 



Table 5 - continued 



Transition 



A (A) Ir/I 



I/I{Uf3) 









(1) 


(2) 


(3) 


(4) 


9-^3 3po 




1 nfin 9 

lUUU.Z 


u.ouz 


n 'is 

U.Oo 






9^3 3po 


-4p 5i 


iUDU.Z 


U.ODZ 


n '^7 
u.o / 






9^,3 3po 
Zp 


4t-» 3c, 
-4p Di 


iUDU.O 


U.ODZ 


n 1 

U. 10 






9r,3 3r~\o_ 
zp 


-•^n 3p^ 
op r2 


1275.0 


0.118 


17.31 






9„3 3r~\o 
zp 


Qr^ 3p- 
op r2 


1 97^ 
IZ / o.o 


nils 

U. J- -LO 


O. lO 






zp Uj^ 


-'^n 3p^ 
op 1^2 


1275.3 


0.118 


0.21 






9^3 3r~)0 
zp Uj^ 


3p^ 


1276 2 


n noA 


'i Q9 

o.yz 






zp 


_Qn 3p, 
op 1 


1276.2 


0.094 


11.62 






zp 


Qt-1 3p„ 
op r^o 


1 97fi s 

IZ / U.O 


n 070 
u.u / y 


(\ 01 
U.U-i- 






zp Ug 


op us 


1343.3 


0.257 


3.30 








op us 


lO^O.U 


n 9^7 


n 'iQ 
u.oy 








op 1-^2 


1345.1 


0.176 


0.68 






9-^3 3pjo 


op JJ2 


1 'iA^ 'i 


n 1 7fi 

U. 1 / u 


n 

U.OO 






9„3 3r)o_ 
zp 


op 1-^2 


1345.3 


0.176 


2.63 






zp Uj^ 


op ui 


1346 4 


n 1 fi7 


1 sn 






zp 


op ui 


1346.4 


0.167 


0.66 






9„3 3po 
zp 


Qr^ 3p- 

op r^2 


1 fi97 'i 
lUZ i .o 


nils 

U. 1 -i-O 


U. -LO 






9t-,3 3po, 
zp ^2 


Qt-> 3p. 
vSp -T 2 


1 fi97 A 
lUZ i 


n 1 1 s 

U. -L J_0 


n A(\ 

U.4:U 






9t-,3 3po 


Qn 3p, 

-op ri 


1 f^oa Q 
iozo.y 


n noA 


n "^n 
u.ou 






9r,3 3po 

zp 


Qri 3p, 

op iri 


1 fi9)R Q 


n OQA 

u.uy^: 


n 07 
u.u / 






9r-,3 3po 

Zp i-Q- 


Qr, 3p, 

-op r^i 


1 ^i9Q 1 
iOZy. i 


n OQA 
u.uy4 


nil 

U. 1 ± 






9r,3 3po_ 
zp r-j^ 


-'in 3p- 
op r-Q 


1629.8 


0.079 


0.18 






n 3 3po 
zp f-^ 


"in 3c, 

-op oi 


1 f\7^ 7 
iU 1 0. 1 


n nsA 

U.Uc54 


^ A'i 
0.40 






9t-.3 3po 
zp ^2 


■iri 3o^ 
op oi 


1 ^7^^ S 
lU I o.o 


n nsA 


S S7 
o.o 1 






9„3 3po 
Zp 


'in 3c, 

-op oi 


1 fi7'^ Q 

ID / o.y 


n nsA 

U.Uo4 


1 S"^ 






9t-,3 3po 


Qr» 3no 
vSp iJ3 


1 7An 'i 


n 9^7 

U.ZO 1 


y.oy 






9f-.3 3po 
Zp f^- 


-op -L'2 


1 7A'i 9 
i / 40. Z 


n 1 7fi 


9 "^Q 

z.oy 






9t-.3 3po_ 
zp 


_Qri 3no 
op iJ2 


1 7A'i 9 

1 / "4:0. Z 


n 1 7fi 

U. -L / U 


7 '^9 
/ .oz 






9t-.3 3po 
zp Jr--]^ 


Qri 3ni 
op JJl 


1 7A^ n 


n 1 fi7 


2 54 






9t-.3 3po_ 
zp ^2 


Qri 3ni 
Op Ul 


1 7A^ 1 

1 / "4:0. 1 


n 1 (\7 

U. -LU 1 


n 1 7 

U. -L i 






9„3 3po 
^P ^0 


op Ul 


1 7A^ 'i 


n 1 fi7 

U. -LU 1 


'i A^ 

0.4:0 






Qo 3po 


An 3c, 
4p 


lOOU.O 


n ■^f^9 
u.ouz 


n 01 

U.U-i- 






Qo 3pO 

Ob 1 ^ 


An 3c^ 

4p oi 


1 ^'i^ fi 

lOO-i-.U 


n '^fi9 
u.ouz 


n 09 

U.UZ 






3po 


An 3c, 
^P '-'1 


1 R'i(\ 9 
lOOU.Z 


u.ouz 


n ofi 

U.UU 






9^3 Ipjo 
zp iJ2 


op 1^2 


oozy. 1 


n A^A 


n 01 

U.U-L 






Qo 3po 


op iJ2 


oyoo.o 


n A^A 

U.'4:0'l: 


n ofi 

U.UU 






Oo i ^ 


op U2 


oyyo.u 


n A^A 


n fin 

U.UU 




1 


Op Jr^i 


Ofi 3po 
0(-l ± 


4114 3 


n ni A 

U.U-l-4: 


n 01 






op rl 


•^rl 3no 
u,2 


4:0 1 O.U 


U.UOO 


n O'i 
u.uo 






op r^i 


Oil -L^^ 


4:0 / y.u 


U.U-i-U 


n 01 

U.U-i- 






'^rx 3-n_ 

op Ul 


-"if] 3po 


AA^iQ Q 

'±4:oy.y 


n 01 A 

U.U-i-^: 


n 1 ^ 

U. -LO 






3p ^Di- 


-3d 


4465.5 


0.014 


0.09 


1.5? 


- 


3p 3D2- 


-3d 


4477.7 


0.014 


0.34 






3p 3D2- 


-3d 


4488.1 


0.033 


0.06 






3p 3D3- 


-3d 


4507.6 


0.033 


0.45 






3p ^Pi- 


-3d 


4564.8 


0.467 


0.02 






3s ^P°- 


-3p ^P2 


4601.5 


0.118 


2.01 


1.5 


1.3 


3s 3pg- 


-3p ^Pi 


4607.2 


0.094 


2.11 


5.7? 


4.2 


3s 3p°- 


-3p 3pi 


4613.9 


0.094 


1.46 


0.9? 


1.0 


3s ^PJ- 


-3p 3Po 


4621.4 


0.079 


2.38 


1.8? 


1.6 


3s 3p°- 


-3p ^P2 


4630.5 


0.118 


6.57 


4.6 


4.8 


3s 3p°- 


-3p ^Pi 


4643.1 


0.094 


2.89 


2.2 


1.5 


3s ip°- 


-3p 3P2 


4654.5 


0.118 


0.21 






3s ip°- 


-3p 3pi 


4667.2 


0.094 


0.19 






3s ipj- 


-3p 3Po 


4674.9 


0.079 


0.26 







_L1 dliOl LrltJli 




T 11 




111 yap) 




3p ^Di-3d 


4774.2 


0.035 


0.16 


- 


- 


3p ^Di-3d 3d° 


4779.7 


0.016 


1.11 


0.8 


1.1 


3p 3D2-3d 3d§ 


4781.2 


0.096 


0.09 


- 


- 


3p 3D2-3d 3d° 


4788.1 


0.035 


1.26 


1.2 


1.4 


3p 3D2-3d ^DJ 


4793.6 


0.016 


0.34 


1.1? 


- 


3p 3D3-3d 


4803.3 


0.096 


1.45 


1.3 


1.9 


3p 3D3-3d 3D§ 


4810.3 


0.035 


0.24 


- 


- 


2p3 iD°-3p iPi 


4895.1 


0.231 


0.07 


1.7? 


- 


3p 3Si-3d 3pg 


4987.4 


0.014 


0.89 


7.6? 


4.6? 


3p 3Si-3d 3P° 


4994.4 


0.014 


2.65 


1.3? 


1.8 


3p 3Di-3d 3F§ 


5001.1 


0.467 


1.10 


- 


1.2 


3p 3D2-3d 3f° 


5001.5 


0.375 


2.30 


- 


1.8 


3s 3p°-3p 3Si 


5002.7 


0.084 


0.53 


- 


- 


3p 3D3-3d 3f° 


5005.2 


0.944 


1.44 


- 


- 


3p 3Si-3d 3p° 


5007.3 


0.033 


3.23 


- 


- 


3s 3p<j>-3p 3Si 


5010.6 


0.084 


1.38 


- 


- 


3p 3D2-3d 3F§ 


5016.4 


0.467 


0.18 


- 


- 


3p 3D3-3d 3f° 


5025.7 


0.375 


0.23 


- 


- 


3s 3p°-3p 


5045.1 


0.084 


2.13 


- 


1.4 


3s iP'j'-3p 


5073.6 


0.084 


0.16 


- 


- 


3p 3Si-3d 3d§ 


5383.7 


0.035 


0.01 


- 


- 


3p 3Si-3d 


5390.7 


0.016 


0.01 


- 


- 


3p 3Po-3d 3P5 


5452.1 


0.014 


0.28 


- 


- 


3p 3Pi-3d 3po 


5454.2 


0.014 


0.36 


- 


- 


3p 3Pi-3d 3po 


5462.6 


0.014 


0.32 


- 


- 


3p 3Pi-3d 3p° 


5478.1 


0.033 


0.18 


- 


- 


3p 3p2-3d 3P5 


5480.1 


0.014 


0.41 


- 


- 


3p 3P2-3d 3p° 


5495.7 


0.033 


0.89 


0.7 


0.5? 


3s 3pj-3p 3D2 


5666.6 


0.176 


4.63 


3.1 


2.9 


3s 3pg-3p ^Di 


5676.0 


0.167 


2.32 


1.2 


1.0? 


3s 3p°-3p 


5679.6 


0.257 


6.24 


4.3 


4.3 


3s Spo-Sp 


5686.2 


0.167 


1.52 


0.8 


0.6? 


3s 3p°-3p 


5710.8 


0.176 


1.52 


0.9 


0.9 


3s 3p°-3p 3Di 


5730.7 


0.167 


0.10 


- 


- 


3s ip°-3p 3D2 


5747.3 


0.176 


0.41 


- 


- 


3s ip'j'-3p ^Di 


5767.5 


0.167 


0.19 


- 


- 


3p 3Po-3d 


5927.8 


0.016 


1.15 


0.7? 


1.0? 


3p 3pi-3d 3d° 


5931.8 


0.035 


1.72 


1.4 


2.0 


3p 3pi-3d 


5940.2 


0.016 


0.80 






3p 3p2-3d 3d° 


5941.7 


0.096 


2.05 


1.2 


1.5 


3p 3p2-3d 3d° 


5952.4 


0.035 


0.51 


0.6? 


1.2? 


3p 3P2-3d 


5960.9 


0.016 


0.05 






3s 3p<j>-3p ipi 


6379.6 


0.231 


0.08 






3s ip°-3p iPi 


6482.0 


0.231 


0.38 






3d 3p°-4p 


6810.0 


0.362 


0.04 




0.3? 


3d 3p°-4p 3Si 


6834.1 


0.362 


0.02 






3d 3pg-4p 3Si 


6847.2 


0.362 


0.01 






4s 3p°-4p 3Si 


14645.6 


0.362 


0.01 







(1) Recombination contribution to total intensity. 

(2) Predicted intensities with an H/3 flux of 2.51 erg cm" 
(3) 



(4) 



BVV 



EPG 



which the LS coupling breaks down. The 3d term is au- 
toionizing, and the most likely mechanism in this case is di- 
electronic recombination. IePGI reported lines at AA 6744.42, 
7535.32 and 9016.42 A that would be produced by highly 
excited d and f states of N^. Their identification as Nil, 
however, is uncertain and their production in our model is 
very unlikely. Many more lines are predicted by our calcula- 
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Figure 3. Comparison of intensities observed bv lBVVl (left) and 
[ePG (right panels) with predicted intensities of fluorescence lines 
from p states (triangles) and d states (circles) for two WMbasic 
stellar atmospheres. Open triangles and circles have uncertain 
/obs (marked with "?" in table 0. The broken horizontal line is 
the value of ii = (Jcal/^obs)- The nebular model was calculated 
with NEBU with BM parameters, except for T^g = 39 kK in the 
lower panels, giving column densities of Af(N+) = 7.3 X 10^^ cm~^ 
and 2.6 x 10^^ cm'^ for the T^s = 37 kK and 39 kK atmospheres 
respectively. Note that the observed intensities here and in tablelsl 
are normalized to I(R/3) = 10*. 



tions, but their i nt ensitie s are lower than the weaker lines re- 
ported bv lBVVl or lEPCI. and their intensities are dominated 
by recombination. Intensities for other lines up to n = 8 and 
I = 1 are available from the authors upon request. Most of 
the observed lines with confident measurements fall within 
0.2 dex of the predicted values in the basic model as shown 
in the upper panels of Fig. |21 



5 VARIATION OF PARAMETERS 
5.1 Stellar temperature 

Fig. shows a comparison of observed and predicted inten- 
sities with two WMbasic atmospheres. The reduction in N"*" 
column density produced by the hotter star reduces the pre- 
dicted intensities in general, but the intensities of the lines 
originating from d states suffer a much greater reduction 
than the ones from p states by an order of magnitude. We 
have traced this effect to an important decrease in the model 
atmosphere flux by a factor of ~ 2 at 529.6 A and a factor of 
- 5 at 533.7 A for T^g > 38 kK as shown in Fig. gl As men- 
tioned in section ITTI 3d states are pumped almost entirely 
by absorptions at those two wavelengths. At the same time 
there is an increase of a factor of ~ 2 in the flux at 508. 7A for 
Teff > 37 kK, which is important in the pumping of 3p states, 
and compensates the decrease in the N"*" column density. 




Figure 4. Surface Eddington flux of two WMbasic atmospheres. 
Arrows show wavelengths of the resonant multiplets at 508.7, 
529.6 and 533.7 A. 



1 1^ 1 , 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 

- \ 

: \ \ [Nil] : 

\ \ 

- ^ ^ - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

'- [OII]^ 
^\ 

\ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


: [oiii]: 


- Nil- 

^ 

I 1 I I 1 1 1 r 1 , r 1 , r , , 1 I , 1 I 1 r 


-i 1 , , I I 1 1 I r 1 1 r r , r 1 , , r , 1 r 


36 37 38 39 40 36 37 38 39 40 



(kK) 



Figure 5. Intensities of some nebular diagnos tic and No permit- 
ted lines normalized to the observed values bv lBVVi as a function 
of stellar Teff. Top left: [Nil] 6584 (solid), [NII]5755 (dashed). Top 
right: [Oil] 3729 (solid), [Oil] 3726 (dashed). Bottom left: [OIII] 
5007 (solid), [OIII] 4363 (dashed). Bottom right: Nil permitted 
lines averaged K = (/cal/^obs) (solid), A 5941.65 A(d upper state, 
short dash) and 4630.54 A (p upp er sta te, long dash). Horizontal 
lines are the values observed by lEPGi The nebular model was 
calculated with CLOUDY with the BM parameters, except for a 
varying T^ff. 
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The difference in the pumping rates of d and p states 
for Teff > 38 kK depends only on the shape of the spectrum 
and the contribution of recombination to the population of 
those states. Thus the disagreement between predicted and 
observed intensities of lines from d states with Toff — 39 kK 
shown in Fig. Elpersists when the N abundance is increased 
to 1 X 10~* or when the recombination contribution to the in- 
tensities is i ncreased with the hyp othesis of ultracold plasma 
proposed bv lTsamis et all i2004r) . A plasma temperature of 
2000 K doubles the predicted intensities of lines from 3p 
states, but lines from 3d states increase their intensities in 
lower proportions because fluorescence is more important in 
the population of those states. 

The other critical parameter in the fluorescence line in- 
tensities is the column density. The N"*" column density de- 
creases with Tcff at a much lower rate for Teff > 38 kK as 
reflected in Fig. |K| The Nil fluorescence lines and the [N II] 
lines decrease little for harder spectra due to that persistent 
N"*" concentration, but their different behaviour at lower TeB 
can be understood in terms of the concentration and the 
escape probability concept. As shown in Fig. |21 the fluores- 
cence N II lines form 50% of their intensity in the inner layers 
of the nebula, much closer to the star than the [N ll] lines, 
which are produced in the outer zone. As Teff decreases, 
the concentration and the optical depth of the resonant 
transitions increase, but the intensities of lines from p and 
d states behave differently as shown in Fig. |5] for the lines 
3s ^P2-3p ^P2A 4630.54 and 3p ^P2-3d ^Dj^A 5941.65. Ab- 
sorption transitions that populate the p states have a much 
lower optical depth than the ones pumping the d states. As 
a result the escape probability decreases more for d states 
than for p states with lower Teff, and the pumping due to 
reabsorption of resonant photons for d states increases. 

NEBU tends to give larger column densities than 
CLOUDY, and thus predicts higher intensities. Predicted line 
intensities by nebu in the UV tend to be 30% more intense 
than those given by CLOUDY because nebu does not consider 
internal dust extinction, but predictions of the two codes are 
within 20% of each other in the optical. 



5.2 Kurucz atmospheres 

CLOUDY contains a g rid of low-resolution Kurucz atmo- 
spheres llKurucz|ll99iri that can be readily used as continua 
in our calculations. A comparison of Fig. |3 and [H] shows that 
the differences between the calculated intensities of p and 
d states with Kurucz atmospheres are much smaller than 
the differences with the WMbasic atmospheres because the 
Kurucz atmospheres do not have the structure of the WM- 
basic atmospheres that causes the different absorption rates 

between p and d states. 

Mode ling of Orion with CLOUDY iSaldwin et al.lll99ll 
Il996l . 1200011 has favored stellar temperatures that are lower 
than current spectro-photometric measurements. Our re- 
sults with the WMbasic atmospheres also favor a lower Tctt, 
but Kurucz atmospheres give a better agreement with ob- 
servations because they are softer than other models and 
give a larger N"*" column. 



2 - 38 kK 



— 38 kK 



I I I I I I I I I I I 



Oa 



40 kK 



_L 



_J_ 



I I I I I I I I I I 



sen 
* • * 

o 



40 kK 



_L 



_J_ 



0.5 0.5 

log /(obs) 

Figure 6. Same as Fig. |21 but for two Kurucz atmospheres. 
The N+ column densities are A''(N+) = 3.5 X 10^^ cm"^ and 
1.13 X lO^'^cm-^ for the T^g = 38 kK and 40 kK atmospheres 
respectively. 





- [Oil] ; 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


[0 111] - 

L /'J 


Nil - 

- . / 

. . 1 " 





12 



13 



12 



13 



Figure 7. Same as Fig. |F] but with T(,fl=37kK and varying </<o. 
The average R = Jcal/^obs (iiot shown) is within 3% of the 
4630.54 A line. The fit of equation 1111 (dots) is normalized to 
the 4630.54 A line at logcf>o = 12.9. 
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5.3 Stellar flux 



ACKNOWLEGMENTS 



Unlike most forbidden and recombination lines, fluorescence 
line intensities are more sensitive to changes in the stel- 
lar flux illuminating the nebula. Their intensities increase 
with (pQ up to 10^^'^ cm~^ s~^ and remain nearly constant 
for higher <j)Q. This behavior can be understood in similar 
terms to the curve of growth of the resonant lines. As the 
intensity of the ionizing flux grows, the N"*" column density 
and the optical depth increase, and the cores of the resonant 
lines become saturated. Eqs. |21 and |U] show that the inten- 
sity of the fluorescence lines is approximately proportional 
to the integral along the line of sight of the pumping rate of 
equation ^ times the density of the absorbing state, UgPgj. 
Changing variable from r to to, eliminating constant quan- 
tities and assuming a constant Doppler width, the intensity 
of a fluorescence line will be 

J,(0) / e-"<= / e-^'>'^^''Uix)dxdro . 



where J,^(0) is the stellar continuum at the illuminated face 
of the cloud, and Tc is the continuum opacity. The integra- 
tion over To can be performed exactly if we assume a mean 
value for e~'^'=. For flxed TbS and v, J^(0) is proportional 
to (po, which in turn is proportional to the H/3 flux. There- 
fore the fluorescence line intensity normalized to H/3 must 
be proportional to 



dx 



The integral is proportional to the curve of growth Wijo). 
Fig. |7| shows that the intensity of the lines follows closely a 
flt of the form 



//7(H/?) oc e 



-2.5xlO"^'']V(H+) 



(11) 



where A''(H^) is the column density in cm ^ and W(tq) 
is the curve of growth of the 2p^ ^P2-4s ^P§A 508.697 A line, 
which pumps most of the 4630.54 line. 



6 CONCLUSIONS 

The intensity of the lines in the N ll spectrum of the Orion 
nebula can be explained by fluorescence of the UV radiation 
of 6^ C Ori in the ionized gas. Recombination of N"*"^ con- 
tributes a minor part of the observed intensities of lines from 
3p and 3d levels connected to the ground state. The effec- 
tive temperature of the star must be below 38000 K in order 
to reproduce the observed line intensities with typical ion- 
ization models that are consistent with the forbidden line 
intensities. An increased N abundance does not allow the 
use of a higher star temperature. The e xistence of interven - 
ing ionized material in the foreground dO'Dell et al.lll993l) 
was not considered in our model and may help increase the 
predicted intensities of the N ll lines. Fluorescence does not 
increase the intensity of the lines from 4f levels, and other 
mechanisms must be proposed to explain their strong inten- 
sities with respect to the collisionally excited and fluores- 
cence lines in the Orion nebula. 



The authors are very grateful to Katia Verner and Gary 
Ferland for valuable advice in running CLOUDY. 
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